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Abstract. We have numerically solved the Grad-Shafranov equation for 
axisymmetric static MHD equilibria of matter confined at the polar cap of 
neutron stars. From the equilibrium solutions we explore the stability of 
O " the accretion mounds using the PLUTO MHD code. We find that pressure 

^-> . driven modes disrupt the equilibria beyond a threshold mound mass, foxvcL- 

^1 1 ing dynamic structures, as matter spreads over the neutron star surface. Our 

results show that local variation of magnetic field will significantly affect 
the shape and nature of the cyclotron features observed in the spectra of 
S^ . High Mass X-ray Binaries. 
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1. Introduction 



Neutron stars in binary systems accrete matter from the companion star, channel- 
t^ ■ ling the matter towards the poles. The accreted matter is confined in a mound by 

the polar magnetic field. Distortions in the local magnetic field due to the pressure 
of the accreted matter can significantly affect the cyclotron resonance scattering fea- 
tures (CRSF) formed there. In the long term such field di s tortions may contribute to 
field burial throu gh diamagnetic screening (Rom anilll990tlMelatos & Phinnevll2001 



Payne & Melatos..2004.) . but the extent of this may be limited by MHD instabilities 



(e.g. Litwin et al. 2001). 
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In this presentation, we first present the solutions of the magnetostatic equations 
describing the accretion mound. We show that accreted matter distort the field lines 
from the unloaded dipolar structure, even at heights much larger than that of the 
mound itself. Next, we perturb the static solutions to study the stability of the system 
and the growth of MHD modes. We find that for mounds above a threshold mass, 
MHD instabilities disrupt the equilibria. We discuss the implications of the local field 
distortions on the cyclotron resonance scattering features (CRSF) emitted from such 
systems and the effect of the instabilities on the long term evolution of the system. 

2. Magnetostatic solutions of accretion mounds 

We consider an accretion mound of polar cap radius Rj, ~ 1 km on a neutron star of 
radius ~ 10 km, mass ~ 1.4 Mq and polar surface field strength ~ 10'~ G, typical of 
mounds on HMXB systems. We consider Newtonian gravity of constant acceleration 
g - -gz- We work in a cylindrical coordinate system (r, 9, z) with the origin at the 
magnetic pole and assume axisymmetry around the z axis. By introducing the flux 
function i//(r, z) describing the poloidal flux through a circle of radius r at a given 
height, one c an recast the static Euler equation into the Grad-Shafranov (hereafter 



GS) equation (jMukherjee & Bhattacharvall2012l) 
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Previously, the approximate equation of state for a non-relativistic degenerate Fermi 
gas (with p oc p^^^) was used to solve the GS equation, which is insufficient to describe 
the plasma for large densities near the base (> 10'' g cm"^). In our current work, 

we have used an equation of state: p - (87r/15)meC^ {^) ^f/ ((1 + 16/25x^)'^^j, 
which closely approximates the T - K Fermi plasma (with errors less than ~ 1.5%; 

IJ 1 I 3 \'^"' 

Paczvnskil (ll983l) V Here ^f = ■;^\ ^^'^ „ 1 p'^^ is the Fermi momentum. The dens- 
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ity can be derived from the expression for Fermi momentum obtained after separation 
of variables: 

1/2 

16 yUpOTn 

A\ ^0/Q ;^=-^(Zo(<A)-z) + l (2) 

4 1 32/9 I 15 m^c-^ 

The shape of the mound is specified by a mound height function Zo(i/') marking the 
top of the mound as a function of flux. Solutions thus obtained show large devi- 
ation from dipolar field configuration, even at heights several hundred metres above 
the mound (s ee Fig.[T). CRSF emitted from su ch columns will have complex shapes 
and features (Mukherjee & Bhattacharval l201 2 |). Some sources like V0332+53 sho\y 
broader CRSF with decrease in luminosity (ITsvgankov. Lutovinov & Serben 120101) 



which can be interpreted as the characteristic emission region coming closer to the 
mound where field distortion is larger Current and future X-ray missions with im- 
proved spectral sensitivity like NuSTAR, ASTROSAT, ASTRO-H, LOFT etc wiU be 
crucial to probe the conditions inside such accretion columns. 



MHD of accretion mounds 
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Figure 1. Left: The field lines inside the column for a GS solution of a 55m mound (solid 
lines), compared with undistorted dipolar field (dashed). The dotted line represents the top of 
the mound. Substantial deviation from dipolar fields extends far above the mound surface, up 
to a height ~ 1km. Right: Ratio of strength of the local field to that of an undistorted dipole, 
as a function of height, at different radial distances from the magnetic axis. For r > 700m 
(where pressure gradients are highest), field lines are pushed outside by confined matter causing 
enhancement of field strength, and a decrease in the inner parts. Even at a height ~ 500m, field 
strength differs by more than 10% of dipole value. 

3. MHD instabilities in accretion mounds 



To investigate the presence of MHD instabilities, we perturb the GS equilibrium solu- 
tion and follow the dynamics with the PLUTO MHD code (Mignone et al. 2007). 
Mukherjee et. al 2013 have shown the presence of gravity driven modes through 2D 
axisymmetric simulations. Here we report the results from 3D non-axisymmetric sim- 
ulations of the mounds. We use GS solutions with p oc p^/^ equation of state, for easier 
numerical implementation of the MHD equations. We have performed these simula- 
tions for mounds of different shapes and masses to study the effect on the growth rates 
of the MHD instabilities. Details of the numerical simulations will be presented in a 
forthcoming paper (Mukherjee, Bhattacharya and Mignone, submitted to MNRAS). 

Mounds of larger mass (~ lO^^^M©) with larger field curvatures are highly prone 
to pressure driven instabilities. Finger like channels appear at the radial edges in a few 
Alfven times as matter passes through regions of low magnetic fields (see Fig.|2]). The 
instabilities develop quickly over time scales of milliseconds. However for mounds of 
smaller mass with less field curvature (e.g. 50m mound with mass ~ 6.8 x IO^'-'M©), 
the growth time scales are ten times slower. Eventually we reach a threshold mound 
size (~ 45m mound of mass ~ 5 x 10^'^Mq) which is stable to perturbations. The 
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Figure 2. Left: Cross section of a 70m mound at a height of ~ 30m and t ~ IIa showing the 
density. Ia is the average local Alfven time ~ 2.8 x 10"^^ s. Right: The magnetic field magnitude 
at the same height and time. The finger like channels due to the MHD instabilities are clearly 
seen at the outer radial edges. See the online journal for a colour version of the figure. 



maximum plasma /? (ratio of plasma pressure to magnetic pressure) for a 45m mound 
is ~ 293 which is close to the threshold/? ~ 260 predicted by Litwin et al. (2001) for 
instability. 



Previous solutions of iMelatos & Phinnevi (1200 Ih and iPavne & MelatosI (120041) 
predict large accretion mounds on the neutron stars formed due to continued ac- 
cretion over long time scales. Such large mounds drag the field lines to form local 
screening currents. However MHD instabilities as presented here, will severely limit 
such screening currents from being formed. Presence of instabilities in much smaller 
mound sizes than previous estimates indicates that such MHD processes will play an 
important role in determining the long term evolution of the field and the spread of 
the accreted matter 
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